EPJ manuscript No. 

(will be inserted by the editor) 



A Schelling-class system with spatially-fixed switching agents 
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Abstract We study the behaviour of a Schelling-class system where spatially-fixed switching agents are 
introduced. Agent types are to be interpreted as social groups, and switching as social mobility. Dynamical 
simulations allow us to provide cross-sections of the 3-dimensional phase diagram F{p, t, /), where p is the 
occupation density, / the fraction of agents able to switch from one group to the other, and r is generally 
interpreted as the tolerance of agents toward neighbours of the opposite type. We investigate different 
transitions between segregated and mixed phases both at / = and along lines of increasing /, where 
the nature of the transition changes. We also find that the presence of switching agents in a segregative 
Schelling-type dynamics can induce de-segregation even with values of / lower than i. 
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1 Introduction 

Since Schelling's seminal work ,1; introducing stochastic 
modeling on a checkerboard to gain insights about seg- 
regation phenomena in urban environments, much effort 
has been devoted to develop and understand multi-agent 
systems subject to Schelling dynamics. The existence of 
similarities with certain types of spin-system and liquid- 
solid dynamics [UI31I11[S1[SJ[Z1IH] has been instrumental both 
in drawing theoretical physicists' attention to segregation 
phenomena and in providing the basis for advanced analo- 
gies with existing results and methods from the statisti- 
cal physics toolbox. Many variants of Schelling's initial 
dynamics have been found to exhibit similar properties, 
and the main characteristic, namely the emergence of seg- 
regation, was shown recently to be a universal property 
common to a vast class of systems [SlITU] . 

We investigate in this article the impact of introducing, 
within a Schelling-class dynamics with two types of agents, 
a given fraction of fixed agents able to switch, that is: to 
change spontaneously from one of the two types of agents 
to the other, rather than simply bearing the same type 
throughout the simulation. We depart here from the stan- 
dard context of Schelling-system simulations where agent 
types usually correspond to ethnic groups. Following Clark 
and Fossett's observation that there is now "powerful ev- 
idence that social distance and preference dynamics do 
have the logical capacity to combine with socioeconomic 
inequality between groups to create relatively high levels 
of ethnic and class segregation" [TT] , we wish to focus on 
the latter and regard agent types as social groups defined 
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eg by income levels. Switching from one type to the other 
can therefore be interpreted as social mobility. Our moti- 
vation is twofold: on the theoretical side, this complexifi- 
cation differs from changes in the agents' utility functions 
that have been explored in the literature, and we would 
like to see how it modifies the system's phase diagram. 
On a more heuristic side, the question is whether "social" 
mobility tends to facilitate mixing and de-segregation, and 
the aim is to explore a new variant of residential dynamics, 
with a long-term view to identifying factors that have the 
logical, mathematical capacity to lead to desegregation or 
lesser segregation. 

We start from a most simple instance of a Schelling- 
class system, which meets the general criteria set out in |5] . 
A noticeable feature of the Schelling dynamics in our sys- 
tem is that the moves are "blind" ones, in the sense that 
when given the opportunity to move, an agent, if they do 
indeed move, will go to a site picked uniformly at random 
among vacant sites. Agents are therefore not satisficers, 
let alone maximizers — but this is not a pre-requisite for 
a system to belong to the general class defined in |2] , as we 
understand it. The specifics of the algorithm implement- 
ing our system are detailed in section [51 Let us emphasize 
that the dynamics studied here is a priori different from 
Glauber and Kawasaki dynamics as discussed in [3] , since 
we impose a non-zero vacancy density and do not allow for 
direct exchange between agents. In section [31 we comment 
on the general behaviour of the system and then derive in 
section[31 from numerical simulations, cross-sections of the 
phase diagram. We compare these to the phase diagram 
in [7], before focusing on the transitions that occur when 
the fraction of agents able to switch increases. 
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2 System and notations 

Ours is a square-lattice system, featuring a finite number 
width X height of sites. At any instant, each site is in one 
of four states: 

1. vacant, 

2. occupied by an agent of (pure) type A, 

3. occupied by an agent of (pure) type B, 

4. occupied by an agent of type C, 

where A and B are the main two types, and C corresponds 
to agents that have the ability to display either type A or 
type B. There is a fixed number of agents of each three 
types (these numbers are determined by the occupation 
density p and the fraction of C-agents /). 

We use periodic boundary conditions and discrete-time, 
ordered, asynchronous updating [T^]. Agents are therefore 
examined in turn ; if they are of types A and B, they fol- 
low a Schelling dynamics, viz. they move to a vacant site 
(chosen uniformly at random) with probability ph if they 
are satisfied with their current position and with proba- 
bility Pu if they are not. The satisfaction of an agent of 
type A or _B is either or 1, depending on the proportion 
of their direct neighbours being of the opposite type: if 
this proportion is smaller than a given parameter r (usu- 
ally interpreted as the degree of "tolerance" of the agents 
[1]) the agent is satisfied, otherwise they are not. Agents 
of type C display at each instant one of the two types A 
or B, but do not move in the Schelling dynamics: when 
examined, they switch their displayed type with probabil- 
ity Ps, regardless of their neighbours' types. The number 
of C-type agents is expressed as a fraction / of the total 
number of agents present in the system. Note that simu- 
lation code and results are available onlin^H. 

Table [I] summarizes the simulation parameters, with 
the typical values which they shall be given hereafter, un- 
less otherwise stated, p, r, /, ph, Pu and Ps were defined 
above. Tm is the number of time steps in a single real- 
ization (one time step corresponds to a complete round 
over all agents), and Ne is the number of ensemble re- 
alizations, width and height are the dimensions in num- 
ber of sites of the simulation grid. The values of p and r 
are chosen so that the pure Schelling dynamics leads to 
a segregated state. (Note that in the basic Schelling-type 
dynamics we are using here, only direct neighbours are 
considered, which implies that values of r tend to be dis- 
tinguished only relative to the discrete values: j, ^ and |; 
however, since vacant sites are not counted as neighbours, 
other values of r can occasionnally be distinguished.) The 
probability of moving for an agent unhappy with their 
current neighbourhood is set to a value significantly lower 
than 1 to reflect difficulties {eg economic factors) that 
may prevent unsatisfied agents from moving. The rate 
Ps of switching for C-agents is then chosen to be about 
one order of magnitude smaller than pu to implement dis- 
tinct timescales for the two types of dynamics (of course 
a reverse order of these timescales should also be studied 
later). Finally, the fact that agents who are statisfied with 



Table 1. Simulation parameters 



Parameter 



Value 



P 

T 
f 

Tm 
Ne 

width, height 

Pu 

Ph 

Ps 



G [0.9, 1] 
0.3 

G [0, 1] 

500 

50 

30 

0.2 

10"* 

0.05 



tolerance= 0.30 density=0.90 f=0.20 



no switching 
switcliing starts at t=0 
switcliing starts at t=500 




200 300 
steps 



500 



|https : / /sourcef orge .net /p/phase-py7] 



Figure 1. Ensemble average {x{t)) of the interface density in 
three simulation modes. 



their current location are able to move {ph ^ 0), even at 
a very small rate, suppresses certain simulation artifacts 
(see eg [1]) and can also be interpreted as reflecting the 
possibily that agents make a "wrong" move due to misin- 
formation or, again, economic factors. 

As recalled in [HITlin], Schelling-class systems are ex- 
pected to reach some kind of equilibrium, or at least a 
steady state as far as segregation is concerned. That trans- 
lates into the reaching of an almost constant limit value 
for a variable quantifying segregation within the system. 
Various choices of variables are possible. We shall follow [9] 
and work with the interface density x{t), defined as the 
average over all occupied sites at time t, of the local in- 
terface density, which is itself equal to the ratio of the 
number of neighbours of the opposite type over the total 
number of neighbours (vacant neighbouring sites are not 
counted). The interface density is normalized by one half, 
so that when the two populations A and B are well mixed 
in the system, x{t) is close to 1, whereas it is close to 
when the system exhibits segregation. We shall write (. . .) 
for averages over all ensemble realizations. 
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3 General results 

Let us first examine the influence of a non-zero fraction of 
switcliing agents on the limit value Xoo = VmYt^oo{x{t)) of 
the interface density (time-averaged on the time-steps af- 
ter a steady state has been reached) , with three simulation 
modes: 

1. no switching agents are present; 

2. a given fraction / = 0.2 of switching agents partici- 
pate, being "activated" only after a fixed delay (500 
time steps); 

3. a given fraction / = 0.2 of switching agents partic- 
ipate, and are "activated" at the beginnning of the 
simulation. 

The first case is clear: this is a standard instance of a 
Schelling system. For the second case, note that C-type 
agents may be "de-activated", i.e. ps may be set to so 
that C-agents keep the same displayed type throughout 
the simulation. They thus act as fixed, non-moving agents 
of one or the other of the two types A and B. There- 
fore, even though they are not themselves moving, they 
do participate in the dynamics by influencing the satisfac- 
tion of their neighbours, and by "persisting" in their res- 
idential choice whatever the current composition of their 
neighbourhood. The third case simply corresponds to an 
activation of C-type agents from the beginning of the sim- 
ulation. 

The occupation density is set to 0.9, the tolerance r 
to 0.3. Starting from a well-mixed state (drawn from a 
uniform distribution), one observes that: 

1. in the absence of switching agents, the system quickly 
relaxes to a segregated state with a value of Xoo very 
close to 0, as expected for the chosen values of p and r 
(knowing the phase diagram of a typical Schelling sys- 
tem as in [J); 

2. in the presence of switching agents that start as fixed A 
or B agents for 500 time steps, relaxation starts toward 
a state of interface density much lower than the initial 
one; however, once activation takes place aX t = 500, 
the ensemble average {x{t)) makes a sudden jump to 
reach a higher limit value close to 0.35; 

3. in the full dynamics, i.e. with C-type agents activated 
from t = 0, the system relaxes toward a moderate value 
of Xoo close to 0.35. 

The introduction of switching agents leads to less- 
segregated states, as could be intuitively expected. It will 
be of interest in future work to examine in more details the 
behaviour of the system in the second variant without acti- 
vating the C-agents: does it have a distinct limit interface 
density Xoo compared to the pure Schelling variant, or is 
the relaxation simply going to take much longer? In other 
words, one may wonder whether the Schelling dynamics is 
just slowed down in its segregated effect, or whether the 
presence of fixed agents, even in proportions less than 20 
percent, leads to some noticeable de-segregation. 

The activation of switching, at time t = 500 in the 
second variant and time t = in the third one, entices 



a significant de-segregation phenomenon, leading quickly, 
in both cases, to a limit value Xoo S [0.3,0.4]. As can be 
seen in the middle part of Fig. [2] (a), for r ranging from 
0.25 to 0.75, this limit value of Xoo remains more or less 
the same. Given that the numbers of A and B agents are 
the same, such a value for the interface density could cor- 
respond to a mixture of phases, as described in 1131: eg 
two "pure" clusters separated by a well-mixed area. That 
could be a way of "optimizing" the average neighbourhood 
(but not necessarily the typical neighbourhood) in a resi- 
dential system at fixed r < ^ : indeed, a notable feature of 
Schelling systems is that, when they lie in the segregated 
phase, their average equilibrium interface density is signifi- 
cantly less than the agents' tolerance. This is precisely one 
of the features that made them famous, since one could 
naively expect that when agents "tolerate" up to one half 
of neighbours to be of the type opposite to theirs, the sys- 
tem would in the long run reach some kind of steady state 
in which each agent had a neighbourhood comprising of 
one half of agents from the other type. But that was firstly 
assuming that average neighbourhood and typical neigh- 
bourhood would coincide, and secondly overlooking the 
possibility of a phase transition occuring at some value 
of the tolerance, especially a value smaller than one half 
— which is the case [TI IT^ITS] . 



4 Phase diagram and transitions 

We wish here to study the transitions occuring in the sys- 
tem both when / = and p, t vary (where we should 
retrieve at least part of the results described in [7]), and 
when / increases. Of course, our main interest lies in the 
latter case, and the relative bareness of our Schelling dy- 
namics (with only direct neighbours taken into account 
and moves being blind) is not suited for an in-depth study 
of the system's phase diagram. 

To facilitate further analysis of the system's behaviour 
and identify the transitions, we shall follow [7] and in- 
troduce, on top of the ensemble-averaged interface den- 
sity {x(t)) and its limit Xoo, a quantity analogous to the 
susceptibility of thermodynamical systems, and another 
analogous to the specific heat. The former is defined as 

^ {[x{t)Y) - {x{t))^ ^ 

T ' 

and the latter as 

Cit)^{[Eit)r)~{Eit))\ 

where E is defined by analogy with the energy in the 
Blume- Emery-Griffiths spin model 

E=-Y^c,c,-{2r-l)Y,^.cl 

with the sums being over pairs of neighbours, and Ci = 1 
for a site occupied by an agent of type A, —1 for an agent 
of type B and for a vacant site. We shall write Xoo 
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(a) Interface density 
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(b) Analogue of suceptibility in log scale 
Figure 2. Interface density Xao and analogue of susceptibility Xoo, in terms of p and r as / increases. 
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(a) Interface density 
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(b) Analogue of suceptibility in log scale 

Figure 3. Profile of the interface density Xco (a) and the analogue of the susceptibility Xoo (b) as r varies, for a given fixed 
occupation density p = 0.95. 



and Coo for the time-averaged steady-state values of x(t) 
and C{t). 

Gauvin et al. have argued and explored the validity of 
such an analogy with thermodynamical spin systems [71 
[S]. We shall follow their lead in trying to identify tran- 
sitions by looking for rapid changes in the value of x^o 
accompanied by peaks in Xoo and/or Coo- 

4.1 Transitions at / = 

We look here briefly at the (/ ~ 0) -section of the phase di- 
agram, which corresponds to the first picture on the left in 



the top part of Fig. [2] for the (/, p, T)-diagram. Of course, 
our system is a much simplified version of the Schelling 
system studied in [7] , and in particular we expect changes 
mainly to occur at a couple of discrete values of r: i, i 
and |. Fig. [5] shows that, at high occupation densities, 
our simulation results are indeed compatible with phase 
transitions occuring at t = 0.25 and t = 0.75. This is also 
supported by the behaviour of the susceptibility, as seen 
in the first figure on the left in the bottom part of Fig. [2] 
The change in the value of the interface density at 
T — 0.25 seems more abrupt than the one at r = 0.75, as 
indicated also by the profile of Xoo(t) at fixed p = 0.95 
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(leftmost plot in Fig. |3](a)). Both transitions are marked 
by peaks of the susceptibility (leftmost plot in Fig.[3](b)). 

In [7], at low vacancy densities, the authors identi- 
fied a first-order phase transition around Tc = 0.75 and 
a "frozen" state below tj = 0.4. In between these values, 
their system would eventually reach a state of low inter- 
face density, that is, a segregated state. The transition 
from a state of frozen dynamics to the segregated phase 
was further investigated in [T7] and found to correspond 
to a jamming transition that could be reproduced via a 
patch model and described by deterministic equations. 

In our system, the low-vacancy-density, low-tolerance 
state seems to be separated from the segregated state by 
a more abrupt (possibly first-order) transition. That state 
could be understood as follows: at low r (that is essen- 
tially T = as soon as r < 0.25), every agent is unsat- 
isfied, unless they have only neighbours of the same type 
as theirs. Therefore, given the opportunity to move, they 
will, because moving does not depend on their finding a 
more "welcoming" site (contrary to what is the case in , 
leading to the dynamics freezing). Eventually, blind moves 
lead to a reasonably well-mixed state, for reasons that dif- 
fer from those behind the frozen, well-mixed state in [7j 
where the system simply remains in its initial state which 
happens to be drawn from a well-mixed distribution. 




Figure 4. Phase diagram in terms of / and t, at constant p — 
0.95. Interface density Xoo (top), analogue of susceptibility Xoo 
(middle) and analogue of specific heat Coo (bottom). 



The general shape of the phase diagram at / = and 
high occupation density is very much comparable to that 
of a classical Schelling system as in |7j , and we can there- 
fore proceed to our main point in the next paragraph, 
with / ^ 0. 



4.2 Transitions along / 

We examine now cross sections of the phase diagram at 
constant p, letting /, the fraction of switching agents 
present in the system, increase. Intuitively, if / is suffi- 
ciently large, one expects a well-mixed system, as the sit- 
uation then becomes that of fixed agents switching with 
equal probabilities to the A or B types. This is what one 
can observe in the upper part of Fig. [JJ where indeed the 
mixed phase (high interface density) gains more and more 
ground as / increases. As we have seen in sectional when 
/ = two transitions are noticed, one at r = 0.25 and one 
at r = 0.75. Given that at large values of /, only the well- 
mixed phase survives, there should be a curved transition 
line joining these two points on the r-axis, jutting into 
the well-mixed phase and enclosing a segregated phase. 
This is precisely what one observes in the simulations, as 
pictured in the upper part of Fig. |4l 

The nature of the transition occuring across that line, 
at different points, is of particular interest, all the more 
as it varies with /: at small values of /, one has the 
transitions described in the previous subsection. They are 
compatible with thermodynamical transitions, marked by 
peaks in the analogues of susceptibility and specific heat 
(lower parts of Fig. Near higher values of /, the tran- 
sition line is no longer parallel to the / axis, and thus can 



be crossed via an increase in /. For instance, at fixed oc- 
cupation density p = 0.95 and tolerance r = 0.3, Fig. [5] 
shows the profiles of the interface density, the analogue of 
susceptibility, and the analogue of specific heat as the frac- 
tion of switching agents / is increased from to 1. The 
change in the interface density is not abrupt but peaks 
can be observed in Xoo and Coo around / « 0.25, with a 
finite-size effect [T51[TF| shown in Fig. (5]. 

Further right in the phase diagram shown in the top 
part of Fig. SI near / « 0.5, where the line separating 
segregated and mixed phases curves to join its upper and 
lower parts, a more gradual change from a low to a high 
interface density is observed at intermediate values of the 
tolerance t when / is increased. No peaks are to be noted 
here in the susceptibility or specific heat analogues. This 
seems to indicate that the nature of the transition has 
changed along the line, possibly from thermodynamical 
to simply a smooth, gradual change. 

Coming back to the plot of x^d against / in Fig. [SJ an 
interesting feature lies in the fact that the interface density 
seems to exhibit a maximum. This would mean that a bet- 
ter mixing (as measured by the interface density) could be 
achieved by tuning / to some intermediate value, rather 
than simply setting it equal to 1. In other words, a combi- 
nation of Schelling dynamics (with segregative effects) and 
switching dynamics could bring more heterogeneity in the 
system than pure switching dynamics. However, this could 
be a bias in what the interface density actually measures 
or a bias linked to our sample of initial distributions: when 
/ = 1, where one would naively expect maximal mixing, 
the value of Xoo is liable to any pattern in the distribution 
of vacant sites since every agent remains fixed. Should it 
turn out to be an artefact linked to the initial distribu- 
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Figure 5. Transition along /, at constant p = 0.95 and 
r = 0.3. Interface density x^o and the analogues of suscep- 
tibility Xoo and specific heat Coo are plotted as the fraction / 
of switching agents increases. 
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Figure 6. Finite-size effect. The analogue of specific heat Coo, 
normalized by its initial value, is plotted for various grid sizes, 
at p = 0.95 and r = 0.3 (with an ensemble of realizations other 
than the one used for Fig.[S|. The larger the system, the more 
marked is the peak, as one would expect in thermodynamical 
phase transitions. 

tion, confirmation of the existence of a maximum at some 
intermediate value of / could still be interesting: one does 
have biased initial states in real cities. 



5 Conclusion 

We have examined the numerical behaviour of a system 
driven by a mixture of two dynamics: (i) a standard 
Schelling preference dynamics where agents of two types 
move stochastically across a grid according to their pre- 
ferred composition of neighbourhood — which can be in- 
terpreted as residential choice, and (ii) a dynamics where 
some agents stay fixed but have the ability to switch from 
one type to the other, which can be interpreted as "social" 
mobility. 



Our findings include the observation of a number of 
transitions as the control parameters p (occupation den- 
sity), r (tolerance for neighbours of the opposite type) and 
/ (fraction of agents able to switch) are varied across their 
domains. Elements indicate in some cases thermodynami- 
cal transitions of the first order, in other cases continuous 
transitions, and in some instances dynamical transitions 
or simply gradual change. We will seek to confirm these 
in future work, paying particular attention to the pres- 
ence of critical points. Finer numerical simulations will be 
conducted to this end, especially with a larger neighbour- 
hood for each agent. Whether the system exhibits eventu- 
ally only two large domains in the segregated phase (with 
a correlation length diverging to infinity) or whether it 
lacks the property to do so, should be examined [TSlli]. 
as well as the nature of the final state in the various in- 
stances of the system: is that an equilibrium state or a 
non equilibrium steady-state? Analytical descriptions of 
the system PUIITU] . at least in some cases, should help to 
better identify order and control parameters, and we will 
try to find mean field and/or simple analytical models that 
capture the dynamical onset of (de)segregation phenom- 
ena, in the manner of what was done in [5] for standard 
Schelling systems. 

Another noticeable result of the simulations conducted 
for this article is the fact that, for some given, segregative 
values of the tolerance r and occupation density p, the 
interface density could show a maximum along a line of 
increasing fraction / of switching agents. In other terms, 
it seems that, without intervening on the occupation den- 
sity nor on the agents' tolerance, mixed or at least lesser- 
segregated patterns of occupation can be achieved in the 
presence of a certain fraction of spatially fixed, but socially 
mobile, agents, without the need for that fraction to be 
close to 1. This property will be further investigated, in 
particular the existence, under specific conditions, of such 
a phenomenon at small values of /, leading to the possibil- 
ity of treating the presence of switching agents as noise in 
the system and observing noise- induced transitions pi] . 

Variants of the system introduced here will also be 
envisaged, in particular: 

— asymmetrizing what we have called "social mobility" , 
eg making it much more likely for C-agents displaying 
type B to switch to type A rather than for C-agents 
displaying type A to switch to type B, as might be the 
case in real socioeconomic contexts; 

— allowing C-type agents to move on the grid, according 
to their displayed type, and again possibly in an asym- 
metric manner, eg by allowing C-agents to move when 
displaying type A but not when displaying type B; 

— emphasizing distinct time-scales for social and spatial 
mobilities; 

— refining the Schelling part of the dynamics, with "non- 
blind" moves 14,22,23,231[7]; random (a)synchronous 
updating; feedback loops (possibly local) acting on the 
tolerance of agents 5j 

— working with more than two displayed types [2511261 

[21]. 
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Exploration of those numerous variants could find useful 
leads in the use of real data p7 l l28 l[24ll29l[T3] and in the 
development of interdisciplinary work with sociologists, 
geographers and urban-system scientists ([SD] rightly crit- 
icizes the lack of interdisciplinary work in the literature 
related to the Schelling model). In particular, it will be 
interesting to look primarily for the existence of a sta- 
ble third phase (other than the segregated and well-mixed 
ones), as emphasized in [IB] . 
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